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Abstract. Chiral order of the Josephson-junction ladder with half a flux quantum per plaquette is stud- 
ied by means of the exact diagonalization method. We consider an extreme quantum limit where each 
superconductor grain (order parameter) is represented by S — 1/2 spin. So far, the semi-classical S — > oo 
case, where each spin reduces to a plane rotator, has been considered extensively. We found that in the 
case of S = 1/2, owing to the strong quantum fluctuations, the chiral (vortex lattice) order becomes dis- 
solved except in a region, where attractive intrachain and, to our surprise, repulsive interchain interactions 
both exist. On the contrary, for considerably wide range of parameters, the superconductor (XY) order 
is kept critical. The present results are regarded as a demonstration of the critical phase accompanying 
chiral-symmetry breaking predicted for frustrated XXZ chain field-theoretically. 



PACS. 75.10.Jm Quantized spin models 
studies 



85.25.Cp Josephson devices - 75.40.Mg numerical simulation 



1 Introduction 

By means of a field-theoretical technique, Nersesyan, Gogolin 
and Efiler claimed [fij that the ground state of the XXZ 
spin chain with sufficiently strong next-nearest-neighbor 
interaction is in the chiral phase, where the spin-screw chi- 
rality is broken spontaneously. Their proposal is astonish- 
ing, because the XXX counterpart, that has been studied 
very extensively so far [|]]|||,|]||@] > ^ s known to be in the 
dimer phase, where the translational symmetry is broken 
spontaneously so long as the frustration is strong enough. 
New treatments devised specifically for XXZ were re- 
ported dJ^| so as to support the aforementioned predic- 
tion. Meanwhile, as for other frustrated XXZ chain, namely, 
the three-leg ladder with diagonal interchain interaction, 
Azaria et al. argued JIo| , [ll|] that a unique critical phase 
possessing high central charge (c = 2) might be realized in 
the ground state. According to their argument, the chiral- 
symmetry breaking is significant to stabilize such high-c 
criticality. These recent developments tell that the XXZ 
anisotropy with in-plane frustration may give rise to un- 
expected exotic phases. The above-mentioned proposals 
are all based on field-theoretical descriptions. Therefore, 
in order to support these scenarios, numerical simulations 
should be carried out. Kaburagi et al. employed the exact 
diagonalization and the density-matrix renormalization- 
group methods for the XX model with the next-nearest- 
neighbor coupling |l2|,^3). They concluded that the chi- 
rality is not broken for the S — 1/2 chain, whereas it 
would be broken for the 5 = 1 chain; namely, numerical- 
simulation result for 5=1/2 appears to contradict the 



field-theoretical description. Nevertheless, numerical sim- 
ulation of such frustrated chain itself is a matter of serious 
methodological concern || ; main obstacles are due to the 
exponentially small energy gap and incommensurability 
of the spin-correlation function, which often miss-matches 
the total system sizes. 

In order to avoid such complications that would arise 
from incommensurability, we have investigated the two-leg 
Josephson-junction ladder, 
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Here, an and are the hard-core boson annihilation and 
creation operators, respectively, at the i-th site along the 
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l-th leg. The parameters ty, t±, V\\ and V± control the in- 
trachain hopping amplitude, the interchain hopping am- 
plitude, the nearest-neighbor intrachain attraction and the 
interchain attraction, respectively. The number of bosons 
is set to N = L (half filled). Throughout this paper, we 
choose in as the unit of energy; namely, 



in = 1. 



(3) 



We have imposed the periodic-boundary condition along 
the ladder; ai y L+i = a/,i- The angle <P denotes the gauge 
twist around each plaquette. Because the angle is set to 
<P = 7T, the ladder is subjected to a uniform magnetic field 
of half a flux quantum per plaquette. That is, suppose 
that the bosons are in the superconducting state (this is 
a subtle issue in one dimension), the magnetic flux would 
possibly be quantized so as to form rigid vortex-lattice 
order along the ladder. This order is confirmed to de- 
velop in the preceding semi-classical analysis That 
is, the chiral order of our model has two-unit-cell period- 
icity. For the purpose of studying the stability of chirality 
against the quantum fluctuation numerically, our model is 
far more advantageous than the frustrated spin chain, be- 
cause the spin-screw chiral order of the latter model has 
long-wavelength incommensurate structure, and thereby 
the numerical data suffer from insystematic finite-size- 
scaling behavior. Tendency to the formation of such vortex 
structure may become more transparent, if we transform 
the boson Hamiltonian ([!]) in terms of spin: Through uti- 
lizing the mapping relations between the hard-core boson 
and the spin-1/2 operators, namely, an 



4 = ^1 



and a} u an — 1/2 = S^, the above boson model is mapped 
to the XXZ spin ladder, 
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apart from a constant term. Note that the signs of the 
Xy-component magnetic interactions are different from 
one leg and the other, and thus an in-plane frustration 
does exist. Therefore, vortex-antivortex alignment gets fa- 
vored. The semi-classical version S ~ * oo of the model 
(|) has been studied extensively In the 

limit, spins are represented by rotators, and through path- 
integral mapping, the system reduces to a two-dimensional 
classical rotator model. According to numerical-simulation 
studies jl4 15 T7|, that rotator model appears to be in the 
KT critical phase. To our surprise, surrounded by such 
critical background of in-plane spin components, rigid long- 
range chiral order develops. And so, it is our motivation 



to investigate the stability of the chiral order against the 
quantum fluctuations of S = 1/2. We show in this paper, 
that the chiral order fails to develop except in a limited 
condition, where an attractive intrachain coupling and a 
repulsive interchain coupling are both turned on. This new 
phase is a demonstration of critical phase accompany- 
ing chiral-symmetry breaking predicted field-theoretically 



Let us mention some remarks concerning the present 
ladder model described by either boson representation eq. 
([!]) or spin representation eq. One might wonder that 
the hard-core condition, in other words, the S = 1/2 con- 
dition, would be too restrictive, and the model misses mi- 
croscopic physical ingredients such as the intra-and-inter- 
grain charge capacitances. Actually, setup of our model 
might be rather phenomenological. We stress, however, 
that the one-dimensional XXZ chains belong to the Tomonaga- 
Luttinger- liquid universality class irrespective of S [fLof . 
(Heuristically, it is known that one-dimensional quantum 
systems possessing the U(l) symmetry are flown to that 
universality in many cases.) As a matter of fact, according 
to the proposal |19], the transverse-correlation exponent 
r)± is governed by the compact formula, 



J z 



2irS 



(5) 



(In our notation, J z is given by J z = — V\\.) That is, 
through varying J z , one can cover all possible range of 
rj± with 5=1/2 fixed. Lastly, we mention the possibil- 
ity of the so-called 'spin liquid' state that often arises in 
quantum spin systems; here, the term 'spin liquid' de- 
notes the state with exponentially decaying short-range 
correlation function, we use the term spin liquid that The 
conventional two-leg ladder |^0|,^l) (with non-frustrated 
interchain coupling) is a prototypical system exhibiting 
spin-liquid state. We will show that our frustrated lad- 
der does exhibit spin-liquid phase for a certain parameter 
condition beside the chiral phase. 

This paper is organized as follows. In the next section, 
we explore the ladder model by means of the exact diago- 
nalization method. The last section is devoted to summary 
and discussions. 



2 Numerical results 

In this section, we present numerical results. We carried 
out exact-diagonalization calculations for the system (^) 
with up to TV = 2L = 32 spins. The data are analyzed in 
terms of the finite-size-scaling theory. 

In Fig. 0, we have plotted the Binder parameter J22|,^3| , 
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for the chiral order, 



M = M chiral = ^(-l) 1 [Su x S 2l ] 2 
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for the system (|j) with tj_ = 0.5, Vji = 0.6 and V± var- 
ied. Here, (■ • ■) denotes the ground-state average. Note 
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Fig. 1. Binder parameter for the chiral order (j7J) is plotted for 
tii = 1, t± = 0.5, Vji = 0.6 and various V±. The symbols o, □, 
o, A and < denote the data for the system sizes L — 8, 10, 12, 
14 and 16, respectively. We see that the chiral order develops 
in the repulsive interchain coupling —1.8 < V± < 0. 



that in the boson language, Af c hi ra i measures the stag- 
gered boson current through the rungs. Therefore, it de- 
tects the formation of the vortex-lattice order. Finite-size 
scaling behavior of the Binder parameter contains the in- 
formation whether the order M develops or not ]22] , |23| ] : 
If the order is long (short) ranged, the Binder parameter 
grows (becomes suppressed) through enlarging the sys- 
tem sizes. At critical point^the Binder parameter remains 
scale-invariant. From Fig. yj, we found that the chiral or- 
der develops for the parameter range —1.8 < V± < 0, 
while it is short ranged otherwise. As would become more 
apparent in the succeeding analyses, the strong-coupling 
regions of V± > 2.5 and V± < —1.8 belong to insulator 
phases of different characters, and are thus rather out of 
present concern. For V± > 2.5, for instance, the bosons are 
so cohesive that they may constitute an island in the sea 
of vacant sites; namely, the system becomes 'phase sepa- 
rated.' On the other hand, for V± < —1.8, owing to the 
strong repulsion, a Mott gap opens between the bonding 
and anti-bonding excitation branches. 

Our result shows that the region of the chiral phase is 
limited. Hence, we found that quantum fluctuations are 
dominating so that unlike the semi-classical case, assis- 
tance of the many-body correlations of Vji and V± is vital 



for stabilizing the chiral order. After scanning the parame- 
ter space, we found that the optimal condition lies around 
t± = 0.5, Vji = 0.6 and V± = —1. It is fairly reasonable 
that the optimal intrachain interaction is attractive, be- 
cause attractive interaction enhances the tendency toward 
superconductivity (XY order); see eq. (||), for instance. 
On the contrary, the fact that the repulsive interchain in- 
teraction is favorable sounds astonishing. It is expected 
that the repulsive interchain coupling might enhance the 
particle exchange between the chains. Hence, we see that 
such particle exchange is significant in order to confine 
(pin) the magnetic fluxes at each plaquette. 

In order to confirm the above phase diagram, we have 
calculated the scaled domain- wall energy, 



LAE(L) = L E S (L) 



LE g (L + l) 
L+l 



(8) 



E g (L) denotes the ground-state energy of the system with 
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Fig. 2. Scaled kink energy (^) is plotted for the same param- 
eter range as that of Fig. [j]. The symbols o, □, o, A and < 
denote the data for the system sizes L = 7, 9, 11, 13 and 15, 
respectively. In the region —1.8 < Vx 0, where the chiral 
phase is realized (see Fig. |l|) , the scaled kink energy grows ac- 
tually. On the contrary, in the attractive region V± > 0, the 
scaled kink energy stays scale-invariant. This indicates that a 
critical phase, probably of the XY mode, is realized; see Fig. 



size L, which is supposed to be odd integer. For the sys- 
tem with odd number of plaquettes, one domain wall is 
created in the alternating alignment of vortex-antivortex 
structure. As would be apparent from the definition (||), 
AE measures the extra domain-wall energy cost due to 
the kink. Hence, for ordered phase, the scaled kink en- 
ergy should increase linearly with respect to L, while for 
disordered phase, LAE vanishes exponentially; namely, 
~ e~ L /£ with correlation length £. At critical point, LAE 
should be scale- invariant, because AE is of scaling dimen- 
sion l/L; note the relations AE ~ l/£ ~ 1/L. In Fig. 0, 
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we plotted LAE for the same parameter range as that of 
Fig. |l|. In fact, in the chiral phase —1.8 < V± < 0, which is 
estimated from the above Binder-parameter analysis, we 
observe clear signature of the chiral-domain-wall energy 
cost. On the other hand, in the area < V± < 2.5, LAE 
seems to be scale-invariant. In the former analysis of Fig. 
[|, we have concluded that in the region, the chirality is 
disordered. Hence, we conclude that in < Vj_ < 2.5, 
the XY order, rather than the chiral order, exhibits crit- 
icality. As are mentioned above, the strong-coupling re- 
gions of V± > 2.5 and V± < —1.8 are belonging to in- 
sulator phases with different characters. In terms of the 
pictures presented above, the behaviors of LAE in these 
regions are readily understandable. The blowup of LAE 
for V± > 2.5 is due to a kink (dislocation) formed in the 
island of particles; note that the particle occupation num- 
ber is forced to be half-filled. On the other hand, the rapid 
closure of LAE for V± < —1.8 is precisely due to the fact 
that the system is a Mott insulator. 

Let us turn to the XY order. In Fig. Q we plotted the 
Binder parameter (ra) for the in-plane spontaneous magne- 
tization (XY order), which corresponds to the supercon- 
ductivity (gauge coherence) order parameter in the boson 
language, 

M 2 = M%y = + Spl^ (9) 

lirnj 

for the same parameter range as that of former figures. 
From the plot, we observe that the XY order is kept 



V± < 2.5. We notice that this XY critical phase does 
contain the chiral phase. This feature is quite contrastive 
to that of the semi-classical case S — > oo, where the chi- 
rality is stronger (more stable) than the XY order, and 
thus the chiral phase contains the XY phase |l5| . More- 
over, it should be noted that in the critical region, the 
fixed-point value of the Binder parameter is kept hardly 
changed. In consequence, we found that the XY corre- 
lation function is not influenced very much by Vj_, while 
the chiral sector is affected significantly by V±. These re- 
sults show that those two sectors behave independently; 
in the commonly referred terminology, those two sectors 
are 'separated.' Such the situation where each mode is de- 
scribed by respective low-energy effective theory occurs 
commonly in one-dimensional physics. This point is dis- 
cussed in Section 

Lastly, let us turn our attention to the possibility of 
spin-liquid state. The simplest way to realize the spin- 
liquid state is given just by setting t± to a very large value. 
In that strong-interchain-coupling limit, the bonding and 
anti-bonding excitation branches are separated so that a 
band gap opens and thereby the ground state becomes dis- 
ordered. Other than that rather trivial way, we found that 
the spin-liquid state is accessible from the set of parame- 
ters treated in the former figures |l|j^ just through setting 
Vji repulsive; see Fig. ||, where we plotted the Binder pa- 
rameter for the XY order for t± = 0.5, V± = — 1 and var- 
ious V\\. We see that for attractive coupling < Vj| < 0.7, 
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Fig. 3. Binder parameter for the XY order (j9() is plotted for 
t\\ = 1, tx = 0.5, V\\ = 0.6 and various V±. The symbols o, □, 
o, A and < denote the data for the system sizes L — 8, 10, 12, 
14 and 16, respectively. The XY order appears to be critical 
in the region —1.8 < V± < 2.5. In the region, the fixed point 
value of the Binder parameter seems to be unchanged. Hence, 
we see that the XY order is not influenced by the change of 
V±. 

critical for considerably wide range of parameter — 1.8 < 



Fig. 4. Binder parameter for the XY order (j9|) is plotted for 
tu = 1, tx = 0.5, Vx — —1 and various Vji. The symbols o, 
□ , o, A and < denote the data for the system sizes L = 8, 10, 
12, 14 and 16, respectively. The XY order becomes disordered 
eventually in the repulsive intrachain coupling Vji < and 
driven to spin-liquid phase. 

the XY sector stays critical as is presented before. (The 
rapid suppression of Uxy in V\\ > 0.7 is due to phase sep- 
aration caused by strong attraction.) On the other hand, 
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for repulsive coupling Vj| < 0, the XY order becomes dis- 
ordered eventually and driven to a spin-liquid phase. This 
behavior tells that the interchain coupling is a relevant 
perturbation in. Vji < 0. That feature coincides with that 
of the conventional non-frustrated ladder p^ , ^5| , p6| . We 
discuss this similarity in the next section. 



3 Summary and discussions 

We have diagonalized the Josephson-junction ladder sub- 
jected to the uniform magnetic field of half a flux quan- 
tum per plaquette (0). Unlike the semi-classical case, the 
chiral order suffers from strong disturbances due to quan- 
tum fluctuations. In order to stabilize the chiral order, we 
need to tune carefully the many-body interactions such as 
V]| = 0.6 and Vj_ = — 1. It is surprising that the repulsive 
interchain interaction gives rise to the stabilization of the 
chiral order. This fact indicates that the particle exchange 
across the chains over the rungs, rather than the gauge co- 
herence, is vital for pinning the vortices at each plaquette. 
On the contrary, we found that the XY order (||) is insen- 
sitive to Vj_ ; the XY order is kept critical for considerably 
wide range of parameters. This result indicates that the 
chiral and XY sectors are separated. 

In quantum spin systems, owing to the strong quan- 
tum fluctuations, the so-called spin- liquid state can appear 
pof . The spin-liquid state, in fact, has been the main con- 
cern in the course of studies of the (non-frustrated) two-leg 
ladder pl| . Close to the chiral phase mentioned above, we 
found that a spin-liquid state emerges just through chang- 
ing the sign of Vji . Surprisingly, this behavior coincides 
with that of the conventional ladder; with the bosoniza- 
tion method [^4 25 2(| , it was shown that the interchain 
coupling becomes relevant for Vu < (antiferromagnetic 



intrachain interaction) . Hence, it is suggested that the XY 
mode of our frustrated ladder and that of the conventional 
ladder behave similarly. This similarity might be reason- 
able, if we remember the separation of the chiral and XY 
sectors mentioned above. Hence, it is suggested that the 
present critical phase, as well as that of conventional lad- 
der p^j2^p?| , belongs to the universality class of the cen- 
tral charge c = 1. fj4p5] , p6]] . Direct identification of the 
universality class in terms of the finite-size conformal-field 
theory, for instance, may be exceedingly troublesome, be- 
cause the excitation structure of the XY sector is smeared 
by that of the chiral sector. This would remain for future 
study. 
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